Convective instability induced by nonlocality in nonlinear diffusive systems 
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We consider a large class of nonlinear diffusive systems with nonlocal coupling. By using a non- 
perturbative analytical approach we are able to determine the convective and absolute instabilities of 
all the uniform states of these systems. We find a huge window of convective instability that should 
provide a great opportunity to study experimentally and theoretically noise sustained patterns. 
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There has been a considerable interest in convective 
instabilities and noise induced patterns recently Q. A 
convective instability happens when a state of a nonlin- 
ear system becomes unstable and a localised perturbation 
grows travelling in the system, but eventually decays at 
any point in the laboratory frame. In this regime a state 
different from the original one cannot be established un- 
less it is sustained by noise. Small regions of convective 
instabilities have been predicted and observed in hydro- 
dynamics || , plasma physics Q and optics Q . These are 
due to spatial drift terms modelled by gradients in the 
direct space or by the first few terms of the Taylor expan- 
sion of the dispersion relation in the Fourier space |j] . 

In this letter, instead, we study nonlinear systems 
with nonlocal coupling. Nonlocality is fundamental in 
modelling material response in presence of transport ef- 
fects 0. Even for local material response, nonlocality 
is unavoidable whenever travelling waves emerging from 
a medium arc reflected back to it non-collinearly, there- 
fore coupling any spatial point x with the shifted point 
x + Ax. Indeed this kind of nonlocal coupling is induced 
by any small misalignment in all feedback optical sys- 
tems (see experiments in Refs. 0,0,0)- We consider the 
most general case in which the spatial shift cannot be 
approximated by a gradient (drift) term in a very broad 
class of nonlinear diffusion equations. Equations of this 
type arise in nonlinear systems with diffraction-free op- 
tical feedback Q and are an interesting generalisation of 
more standard nonlinear diffusion equations. We show 
that the stability of all the uniform states of this class of 
systems is determined by a single dispersion relation with 
two parameters. By using a non-perturbative analytical 
approach we are able to analyse this dispersion relation. 
We find a huge convective instability window where noise 
sustained patterns and amplification of perturbations can 
be observed. 

We consider equations of the type 

(d t - %)<f>(x, t) = M^x, t); it) + f 2 {<j>{x + Ax, t); fi) (1) 

where <f> is a real variable, t is in units of the diffusion 
time, x and the spatial shift Aa; are in unit of the diffusion 
length, /i is a control parameter independent on x. f\, /2 
are real functions that can be derived with respect to <j). 
The uniform states 4> m of Eq. JJ) are the solutions of /i + 



fi = and their domains of existence depend upon \i but 
not upon Aa; in the limit of infinitely extended systems. 
The dispersion relation for perturbations exp (cot + ikix) 
of a uniform state of Eq. JJJ is 

; n)e iklAx . 



u = -kj + d,f,fi(4> m ; n) + d rt ,f 2 (4>r, 



(2) 

As a consequence of the term e lk ' Ax , which is present 
in all systems with shift, there are bands of fcj for 
which the real part of the dispersion relation ujr can 
be positive. For d^fi{4>m\ fx) < 0, as in the experi- 
ments in Ref , these bands are within the regions where 
d&fz (0m ; m) cos kiAx > 0. As a result the homogeneous 
solution 4> m is unstable and plane wave perturbations are 
amplified. As the imaginary part of to (phase velocity) is 
in general non null, these waves move across the system. 
A peculiar effect of the non-locality is that for spatially 
localised perturbations the sign of the group velocity and 
of the phase velocity of the most unstable wavenumbcr 
are always opposite. 

When the group velocity is non vanishing, we have 
to consider whether localised perturbations produce ab- 
solute or convective instability, that is, to find whether 

the Green function, f + °° e^^^^'^dki. of the linearised 

n 

equation diverges or vanishes for t — > +oo 8] . To evalu- 
ate the integral, it is convenient to extend analytically 
lo in the complex plane k = kit + iki and apply the 
saddle-point method. An even number (at least four) of 
paths with luj constant start from each saddle point with 
dcu/dk = 0, as shown in Fig.[I] On half of these equiphasc 
paths, the steepest descent, ur decreases fastest; on the 
remaining half, the steepest ascents, wr increases fastest. 
If one can form a closed integration contour with the 
imaginary axis and steepest descents, then the asymp- 
totic value of the integral is given by the values of uj 
at the saddle points on the integration contour ||. This 
method has been extensively used to find the threshold 
between convective and absolute instabilities in systems 
with drift, where the dispersion has in general few sad- 
dle points. However for nonlocal systems the exponential 
term in the dispersion originates always a countable in- 
finity of saddle points. Moreover, the saddles and their 
steepest descents move and can suddenly collide and dis- 
appear as the control parameters change. It is therefore 
essential to study very carefully how the global gcomctri- 
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FIG. 1: Equiphase paths in the complex plane, wi(qn,qi) — 
constant, with w = loAx 2 and q = kAx, for 3 values of 
the control parameter d r f,f2(4>m', (i)Ai 2 = 4 (a), —0.5 (b), 
— 1 (c). The saddles points Si are at the intersections be- 
tween equiphase paths. The symbols show the steepest de- 
scent paths closing the integration contour. 

cal organisation of the saddles and of the equiphase paths 
changes with the control parameters in order to close cor- 
rectly the integration contour. 

In order to simplify the calculations, we define the 
parameters f = 9^/i(0 m ; /x), x = ^/ 2 (^ m ;/x). The 
analytic extension of the dispersion relation is w = 
£Ax 2 + q 2 + \Ax 2 e q , with w — ujAx 2 and q = kAx, 
and we consider only the semi-plane qi > because 
w(q) = w*(q*). Note that the dispersion relation and 
the stability depend only upon the effective parameters 
t;Ax 2 and x' = xAx 2 . For all values of x' there is at 
least a saddle and at most two saddles in the interval 
qi e [0, 7r). We call the leftmost of these saddles so (see 
Fig. QJd). For all n > 0, there are also saddles s„ in the 
n-th interval qj <E (2wr, (2n + l)ir) for x' < (Fig.[I]>c) 
or in the n-th interval qj g ((2n — l)7r, 2ri7r) for \' > 
(Fig. [TJi). In order to close the integration contour en- 
tirely with steepest descents, we need a steepest descent 
path that ends on or is asymptotic to the real axis (qn) 
and another asymptotic to the imaginary axis (qi). These 
steepest descents must either come from the same saddle 
or be connected by other steepest descents. 

From the equation of the equiphase paths that orig- 
inate at the saddle s, 2qjqn + x'e 9H sing/ = Wi(s), we 
can show that the possible asymptotes of the steepest 
descents are the imaginary axis for qi — > +oo and the 
lines qi = (2n + l)ir for x' > an d n > (or the lines 
qi = 2nir for x' < 0) and qn — » +oo. We determine 
the geometrical organisation of the steepest descents by 
using this information, together with the exact determi- 
nation of the interception of the steepest descents with 
the lines qi = nir and with the imaginary axis and the 
fact that, in this problem, two steepest descents with dif- 
ferent phases can intersect only at infinity ^(|- We find 
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FIG. 2: Instabilities diagram: Stable (S), convectively unsta- 
ble (C) and absolutely unstable (A) regions for £ < 0. The 
dashed lines indicate where different bands become unstable. 
The insert shows the well-known diagram [TH obtained when 
instead of a non-local coupling a walk-off term (v) is consid- 
ered. In this case there are no instabilities at all in the lower 
diagram region. 



analytically that the saddle with the smallest phase wi 
(si in Fig.^, and s in Fig.^3 and c) has always the up- 
per steepest descent to the left of all the other steepest 
descents and asymptotic to the imaginary axis. Indeed 
this saddle is always necessary to close the integration 
contour. If so happens to be the saddle with the smallest 
phase, then the steepest descents of so close the inte- 
gration contour (Fig. Q]>c) because so either lies on the 
real axis or has another steepest descent asymptotic to 
the real axis. If instead there is another saddle s ni with 
w;(s ni ) < wi(sq) (for instance the saddle s\ in Fig. ^i), 
than the upper steepest descent from so remains below 
s m and is connected to the steepest descent from s ni (in 
Fig. 2Ji these two steepest descents are asymptotically 
connected at qi = n). In this case we use a steepest de- 
scent from s ni to reach values of qi above s ni itself and 
a steepest descent from so to reach the real axis. If s ni is 
the saddle with the smallest phase, then the steepest de- 
scents of so and s ni close the integration contour (symbol 
lines in Fig.^Ji), otherwise we will have to include also the 
first saddle, s„ 2 , with n-i > n\ and wi(s n2 ) < Wi(s ni ). 
In order to find the saddles to be used to close the inte- 
gration contour, we repeat this process including all the 
saddles whose phase is the minimum of the phases of the 
saddles below them. For each finite x', this procedure 
gives a finite set of saddles P x > 0] . Examples of differ- 
ent integration contour are shown in Fig.Q as symbols 
lines. 

From Eq. J2J we identify the region Wr(0, qi) < in 
which the homogeneous solutions are stable. This is the 
region S in Fig.|2land its contour is the convective thresh- 
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old. If wr(0, qi) > for some qi and if at least one sad- 
dle s n £ P x ' has wr(s„) > 0, the instability is absolute. 
However, we need to find only a small subset of saddles in 
P x > to determine the stability because Wr(s„) > only 
if there is instability with u>i?,(0, qj) > for qi in the n-th 
interval |l3j . Therefore, for each ^, we can determine 
the nature of the instability by finding the values of w 
at the saddles in P x > that are in bands with instability 
threshold x'n below the threshold x'o where wn(so) = 0. 
The absolute threshold %' a is then the value of \' where 



max{w R (s n )\s n g P x >,x"n < x'o) = °- 



(3) 



The importance of the determination of P x ' and of 
Eq. (PJ) is twofold. On one hand they guarantee that we 
can apply the method of steepest descents by properly 
closing the integration contour; on the other hand, they 
allow us to find the absolute threshold simply by inspec- 
tion of w at a finite number of saddles. This is remarkable 
in view of the infinite number of saddles produced by the 
shift. Moreover, P x > is the same for all systems with the 
same x' m the class considered because the global geo- 
metrical organisation of the saddles and of their steepest 
descents does not depend on £. 

Applying this technique we obtain the instabilities di- 
agram in Fig. [21 valid for any £ < (including the typical 
case of a linear damping term £ = —1). For x < the 
instability bands have qj ^ and the lowest convective 
threshold is very far from the absolute threshold. For 
X > 0, the instability with respect to perturbations with 
qi = is absolute for £ = and convective for £ ^ 0, 
with the convective instability windows increasing as £ 
decreases. Comparing the instabilities diagram for dif- 
fusive problems with finite shift and with with walk-off 
(Fig. EJ we note that the whole modulation insta- 
bility region for negative values of x is a specific effect 
of a finite shift, as was already recognised by What 
our analysis reveals by direct calculation of the absolute 
threshold is that in this parameter region the system is 
mainly convectively unstable and shows noise sustained 
modulated patterns. Only for very negative values of 
x/\£,\ the absolute threshold is crossed. 

To show the general applicability of the stability anal- 
ysis presented in this letter we consider in the following 
two examples: the Ginzburg-Landau and the saturable 
nonlinear equations, both with nonlocal nonlinear terms. 
A nonlocal Ginzburg-Landau (NGL) equation is obtained 
by Eq. (QJ when 



(4) 



>± = 



The three uniform states (f>o = (for any /i) and 
±^//j, — 1 (for fj, > 1) are associated to the parameters 
Co = £± = -1 and xo = Mi X± = 3 - %l (Fig. EH)- 
Given the relation between xo,± and /U, Fig. provides 
the linear stability diagram of each state of the NGL 
equation. In particular, x± is always less than 1 so that 
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FIG. 3: a) Xo(aO an d X±(aO f° r the NGL equation, b) Gaus- 
sian perturbation (continuous line) of the vanishing state </>o 
and evolved state for n = 1.03 (dotted line) and for fi = 1.09 
(dashed line). Shift Ax = 0.48. The instability for u = 1.03 is 
convective (local decay, even if the maximum of the perturba- 
tion grows), while for /i = 1.09 it is absolute (the perturbation 
grows also locally). After a longer transient the perturbation 
ends up colliding with the uniform stable states </>+ or 
Simulations with space discretization 0.12, 4096 grid points 
and time step 0.001. 




FIG. 4: a) Xo(m) an d X±(m) f° r the AS equation, b) Noise 
sustained pattern for Xo/|£| — A* = "5 and Aa- = 1.92 in the 
NAS with an additive Gaussian white noise of amplitude 0.01. 
Simulations with Dirichlet boundary conditions and same nu- 
merical discretization of Fig. E| 



only the lower part of the diagram in Fig. El describes <fi± 
instabilities, while for the state <fio = all the instabilities 
shown in the diagram arise by varying and the shift Ax. 
Without shift the state </>q becomes unstable for /i > 
1 where the system evolves to the homogeneous states 
<j> + and 0_ (connected by fronts). With a finite shift for 
fi > 1 a region of convective instability arises before the 
absolute threshold. Simulations of the NGL equation 
with small values of \i show instabilities of the states 4>q : ± 
in good agreement with our theoretical analysis. Fig. EH> 
is an illustrative example of the typical evolution of an 
initial perturbation of the unstable state 4>o, below and 
above the predicted absolute threshold. 

We now consider the nonlocal saturable nonlinear 
(NSN) equation obtained by Eq. when 



fi = 



h = M 



1 



We find the same three uniform states 
— 1 as before, but now xo = Mi X± - 



2-u 



(5) 



(Fig. EH). 
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FIG. 5: Temporal evolution of an initial Gaussian in x — 
perturbing the vanishing state cj>o. a) For fj, = — 14 the 
perturbation evolves toward a modulated wave-packet that 
eventually leaves the system (convective instability), b) For 
(X = — 18 the right front moves to the positive x space region, 
and finally the modulated solution occupies the whole system 
(absolute instability). Simulation of the deterministic NSN 
equation with same parameters as in Fig. UJj. 

Again, knowing Xo,±(A*) an d £ = — 1, Fig. Ogives the 
thresholds of the NSN in terms of the specific parameters 
fi and Ax. We have seen before that for negative values 
of x (here \o = A* < 0) new modulation instabilities are 
predicted (both convective and absolute), existing only 
for non- vanishing shifts and not observed in systems with 
drift. In Fig.^J) we show an example of a noise sustained 
stripe pattern, as predicted by our theoretical analysis. 

The NSN equation has the advantage of not developing 
divergences even for large values of the control parameter 
due to the saturable nonlinearity. Therefore it was possi- 
ble to check the instability diagram even for large nega- 
tive values of xo, decreasing fi from —14 to —18 (Fig. [SJ). 
Due to the modulational character of the instability an 
initial Gaussian perturbation does not evolves in the sim- 
ple way shown in Fig. \Qp. Nevertheless from the tempo- 
ral evolution of the fronts it is still possible to distinguish 
without ambiguity between convective (Fig. |Sp,) and ab- 
solute (Fig.03) instabilities. It is also interesting to note 
that the wave-packet shows a clear asymmetry between 
the left and the right edges, with the wavenumbers with 
lower phase velocity in the leading edge and those with 
higher phase velocity in the trailing edge. This effect 
results from the symmetry breaking caused by the shift 
and the opposite sign of the group velocity and the phase 
velocities of most wavenumbers. 

In conclusion we have identified nonlocality as a new 
mechanism leading to convective instabilities in a large 
class of systems. We have shown how to determine the 



absolute threshold by finding the values of w at a fi- 
nite number of saddles which are selected out of an infi- 
nite number by the procedure described here. Our non 
perturbative approach allows us to analyse situations in 
which there are several bands of unstable wavenumbers 
and the absolute instability is very far from the lowest 
convective threshold. Theoretical predictions have been 
confirmed by numerical simulation of the dynamics of two 
prototype models. We expect that experiments in nonlo- 
cal systems will exhibit a dynamics dominated by noise 
and boundary effects in a very large region of control 
parameters. This opens the possibility of new and excit- 
ing research in the fundamental properties of extended 
nonlinear systems. 
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